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A Novel Shape Parameterization Approach 

Jamshid A. Samareh* 

NASA Langley Research Center, Hampton, VA 23681 

This paper presents a novel parameterization approach for complex shapes suitable 
for a multidisciplinary design optimization application. The approach consists of two 
basic concepts: (1) parameterizing the shape perturbations rather than the geometry 
itself and (2) performing the shape deformation by means of the soft objects animation 
algorithms used in computer graphics. Because the formulation presented in this paper 
is independent of grid topology, we can treat computational fluid dynamics and finite ele- 
ment grids in a similar manner. The proposed approach is simple, compact, and efficient. 

Also, the analytical sensitivity derivatives are easily computed for use in a gradient-based 
optimization. This algorithm is suitable for low-fidelity (e.g., linear aerodynamics and 
equivalent laminated plate structures) and high-fidelity analysis tools (e.g., nonlinear 
computational fluid dynamics and detailed finite element modeling). This paper contains 
the implementation details of parameterizing for planform, twist, dihedral, thickness, and 
camber. The results are presented for a multidisciplinary design optimization application 
consisting of nonlinear computational fluid dynamics, detailed computational structural 
mechanics, performance, and a simple propulsion module. 


Nomenclature 

A wing area 

SR wing aspect ratio 

B Bernstein polynomial 

b wing span 

C chord 

c camber 

cl degree 

e scale factor for twist and shearing 

N B-spline basis function 

n normal vector 

P coordinates of NURBS control point 

R coordinates of deformed model 

r coordinates of baseline model 

S shearing vector 

T twist plane 

t thickness 

u parameter coordinate 

v design variable vector 

W NURBS weights 

X, Y, Z Cartesian coordinates of deformed model 
x, y, z Cartesian coordinates of baseline model 
a angle of attack, deg 

A total deformation 

deformation 

9 twist angle, deg 

A leading edge sweep angle, deg 

A wing taper ratio 

£, r), £ coordinates of deformation object 
p twist radius 

Subscripts 

ca camber 

*Research Scientist, Multidisciplinary Optimization Branch, 
Mail Stop 159, j.a.samareh@larc.nasa.gov. 


/, J, I\ total numbers of control points 
i,j, k indices for NURBS control point 

L wing lower surface 

le leading edge 

m center 

p degree of B-spline basis function in i direction 

pi planform 

q degree of B-spline basis function in ij direction 

r root 

sh shear 

te trailing edge 

t tip 

th thickness 

tw twist 

U wing upper surface 

Superscripts 

T transpose of the matrix 

Introduction 

M ULTIDISCIPLINARY design optimization 
(MDO) methodology seeks to exploit the 
synergism of mutually interacting phenomena to 
create improved designs. An MDO process commonly 
involves sizing, topology and shape design variables. 
Multidisciplinary shape optimization (MSO) finds the 
optimum shape for a given structural layout. It is 
a challenging task to perform MSO for a complete 
airplane configuration with high-fidelity analysis 
tools. The analysis models, also referred to as grids 
or meshes, are based on some or all of the airplane 
components, such as skin, ribs, and spars. The aero- 
dynamic analysis uses the detailed definition of the 
skin, also referred to as the outer mold line (OML), 
whereas the computational structural mechanics 
(CSM) models use all components. Generally, the 
structural model requires a relatively coarse grid, but 
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it must handle very complex internal and external 
geometries. In contrast, the computational fluid 
dynamics (CFD) grid is a very fine one, but it only 
needs to model the external geometry. The MSO of an 
airplane must treat not only the wing skin, fuselage, 
flaps, nacelles, and pylons, but also the internal struc- 
tural elements such as spars and ribs (see Fig. 1). The 
treatment of internal structural elements is especially 
important for detailed finite element (FE) analysis. 
For a high-fidelity MSO process to be successful, the 
process must be based on a compact and effective set 
of design variables that yields a feasible configuration. 
For more details, readers are referred to an overview 
paper by this author 1 on geometry modeling and grid 
generation for design and optimization. 

The model parameterization is the first step for an 
MSO process. Over the past several decades, shape 
optimization has been successfully applied for two- 
dimensional and simple three-dimensional configura- 
tions. The recent advances in computer hardware and 
software have made MSO applications more feasible 
for complex configurations. An important ingredient 
of aerodynamics shape optimization is the availability 
of a model parameterized with respect to the aero- 
dynamic parameters such as planform, twist, shear, 
camber, and thickness. The parameterization tech- 
niques can be divided into the following categories: 2 
discrete, polynomial and spline, computer-aided de- 
sign (CAD), analytical, and deformation. Readers are 
referred to reports by Haftka, 3 Ding, 4 and Samareh 2 
for surveys of shape optimization and parameteriza- 
tion. 

In a multidisciplinary application, the parameteriza- 
tion must be compatible and adaptable to various anal- 
ysis tools ranging from low-fidelity tools, such as linear 
aerodynamics and equivalent laminated plate struc- 
tures, to high-fidelity tools, such as nonlinear CFD and 
detailed CSM codes. Creation of CFD and CSM grids 
is time-consuming and costly for a full airplane model: 
it takes several months to develop detailed CSM and 
CFD grids based on a CAD model. To fit, the MSO 
process into the product development cycle times, the 
MSO must rely on the parameterization of the analy- 
sis grids. For a multidisciplinary problem, the process 
must also use a geometry model and parameteriza- 
tion consistently across all disciplines. For use with 


gradient-based optimization, the geometry model must 
provide accurate sensitivity derivatives of the analysis 
model with respect to design variables. 

This paper presents an approach for shape parame- 
terization suitable for a multidisciplinary design opti- 
mization application. The approach consists of two 
basic concepts. The first concept is based on pa- 
rameterizing the shape perturbation rather than the 
geometry itself. The second concept is based on us- 
ing the soft object animation 5 (SOA) algorithms for 
shape parameterization. The combined algorithm, ini- 
tially introduced by this author, 6 was successfully im- 
plemented for aerodynamic shape optimization with 
analytical sensitivity with structured grid 7,8 and un- 
structured 9 grid CFD codes. 

Parameterizing the Shape 
Perturbations 

At first sight parameterization by splines may seem 
to be a viable approach for shape parameterization. 
The spline representation uses a set of control points 
to define any shape. These control points could be 
used as design variables for optimization. Typically 
over a hundred control points are required to define 
an airfoil section and over 20 airfoil sections to define 
a conventional wing. This requirement results in over 
two thousand control points (i.e., six thousand shape 
design variables) for a simple wing. The number of 
control points is even larger for a complete airplane 
model created with a commercial CAD system. The 
large number of control points is needed more for the 
accuracy than for the complexity. 

Even if we could afford to use a large number of 
design variables, the automatic regeneration of analy- 
sis models (e.g., CSM and CFD grids) is not possible 
with the current technology. For example, it takes 
several months to create an accurate CSM model of 
an airplane. Also, traditional shape parameterization 
processes parameterize only the OML and are inef- 
fective in parameterizing internal components such as 
spars, ribs, stiffeners, and fuel tanks (see Fig. 1). 

It is possible to use any shape (e.g., a sphere) as the 
initial wing definition, allowing the optimizer to find 
the optimum wing shape; however, it is not a com- 
mon practice. Typically, the optimization starts with 
an existing wing design, and the goal is to improve 
the wing performance by using numerical optimiza- 
tion. The geometry changes (perturbations) between 
initial and optimized wing are very small, 10,11 but the 
difference in wing performance can be substantial. An 
effective way to reduce the number of shape design 
variables is to parameterize the shape perturbations 
instead of parameterizing the shape itself. Through- 
out the optimization cycles, the analysis grid can be 
updated as 
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Fig. 2 A typical MSO process. 


R(v) = r + A R(v) (1) 

where r is the baseline grid, R is the deformed (per- 
turbed) grid, A R is the change (perturbation), and 
v is the design variable vector. The change, A.R, is 
a combination of changes in thickness, camber, twist, 
shear, and planform: 

A R = SRth + SR. C a + d'f?tw + (2) 

It takes far fewer design variables to parameterize the 
shape perturbation A R than r itself. 

Figures 2-3 contrast the typical and modified MSO 
processes. In a typical MSO process (Fig. 2), a geom- 
etry modeler perturbs the baseline geometry model. 
Because automatic grid generation tools are not avail- 
able for all disciplines, it would be very difficult to 
automate this MSO process. In contrast, the modi- 
fied MSO process (Fig. 3) relies on parameterizing the 
baseline grids, hence making it possible to automate 
the entire MSO process. 

Soft Object Animation 

The field of SOA in computer graphics 5 provides 
algorithms for morphing images 12 and deforming mod- 
els. 13,14 These algorithms are powerful tools for modi- 
fying shapes: they use a high-level shape deformation, 
as opposed to manipulation of lower level geometric 
entities. Hall 12 presented an algorithm and provided 
computer codes for morphing images. The defor- 
mation algorithms are suitable for deforming models 
represented by either a set of polygons or a set of para- 
metric curves and surfaces. The SOA algorithms treat 
the model as rubber that can be twisted, bent, tapered, 
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Fig. 3 The modified MSO process. 

compressed, or expanded, while retaining its topology. 
This is ideal for parameterizing airplane models that 
have external skin as well as internal components (e.g., 
see Fig. 1). The SOA algorithms relate vertices of an 
analysis model (grid) to a small number of design vari- 
ables. Consequently, the SOA algorithms can serve as 
the basis for an efficient shape parameterization tech- 
nique. 

Barr 13 presented a deformation approach in the con- 
text of physically based modeling. This approach uses 
physical simulation to obtain realistic shape and mo- 
tions and is based on operations such as translation, 
rotation, and scaling. With this algorithm, the defor- 
mation is achieved by moving the vertices of a polygon 
model or the control points of a parametric curve 
and surface. Sederberg and Parry 14 presented another 
approach for deformation based on the free-form de- 
formation (FFD) algorithm that operates on the whole 
space regardless of the representation of the deformed 
objects embedded in the space. The algorithm allows 
a user to manipulate the control points of trivariate 
Bezier volumes. The disadvantage of FFD is that the 
design variables may have no physical significance for 
the design engineers. This drawback makes it difficult 
to select an effective and compact set of design vari- 
ables. This report presents a set of modifications to 
the original SOA algorithms to alleviate this and other 
drawbacks. 

For the modified SOA algorithms presented in the 
next several sections, implementation will include the 
following common set of steps: 

1. Select an appropriate deformation technique and 
object. This defines the forward mapping from 
the deformation object coordinate system (£, i], (j 
to the baseline grid coordinate system ( x , y, z). 

2. Establish a backward mapping from the baseline 
grid coordinate system ( x , y, z) to the deforma- 
tion object coordinate system (£, r], (j. The £, y, C 
mapping parameters are fixed and are indepen- 
dent of the shape perturbations. This is a prepro- 
cessing step that is required only once. 
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3. Perturb the control parameters (design variables) 
defining the deformation object. 

4. Evaluate the grid perturbation (A R) using the 

rj, £ parameters. 

The following sections provide recipes for using SOA 
algorithms for parameterizing airplane models for 
thickness, camber, twist, shear, and planform changes. 

Thickness and Camber 

We use a nonuniform rational B-spline (NURBS) 
representation as the deformation object for thickness 
and camber parameterization. The NURBS represen- 
tation combines the desirable properties of National 
Advisory Committee for Aeronautics (NACA) defini- 
tion 15 and spline techniques, and it does not deterio- 
rate or destroy the smoothness of the initial geometry. 
Readers should consult the textbook by Farin 16 for a 
detailed discussion on NURBS. 

The changes in thickness and camber are repre- 
sented by 

E ;,p(£) E 

SRt*(Z, v) = ^7 ( 3 ) 

E NiAQZNjMWij 

i=0 j =0 

E Ni,p{£) E Nj : q( r })\'Vi'j P c j 

SRc,((, V) = —7 ^7 (4) 

E 

2=0 j =0 

where P t h, , and P c&i ■ are control points (forming a 
control surface) for thickness and camber, Wij are the 
weights, and Ni iP and Nj >q are the p and q degree 
B-spline basis functions defined on the nonperiodic 
and nonuniform knot vectors. Figures 4-5 show the 
NURBS control points in (£, rj) and ( x , y , r) coor- 
dinate systems, respectively. The control points and 
weights could be used as design variables. 

The NURBS representation has several important 
properties for design and optimization. A NURBS 
curve of order p, having no multiple interior knots, 
is p — 2 differentiable. As a result, the NURBS repre- 
sentation can handle a complex deformation and still 
maintain smooth surface curvature. Readers are re- 
ferred to the textbook by Farin 16 for details on the 
properties of NURBS representation. The control 
points are the coefficients of the basis functions, but 
the smoothness is controlled by the basis functions, 
not the control points. The NURBS representation is 
local in nature, allowing the surface to be deformed lo- 
cally, hence leaving the rest of the surface unchanged. 
Equations 3 and 4 serve as the forward mapping be- 
tween the thickness and camber design variables and 
the grid perturbation ( (fifth: SR ca ). 



Fig. 4 Thickness and camber definitions in wing 
coordinate system. 



Fig. 5 Thickness and camber definitions in a;, y , 
and z coordinate system. 


The next step is to establish the backward mapping 
from the deformation object (i.e. , NURBS surface) 
coordinates (£,??) to the baseline model coordinates 
(x,y, z). The wing coordinate system — percent chord 
and span — is a good candidate. The percentage chord, 
%c, is used for and the spanwise location, y. is used 
for rj. 


t = %c, r) = y (5) 

To calculate %c, we need to determine the wing chord 
at each y station. The baseline CAD model provides 
the leading edge (Ri e (r])), trailing edge (R te (r})), wing 
center R m {ri), and normal vector defining the airfoil 
plane T(p) as shown in Fig. 6. The curve defining the 
wing center does not have to be at the center of the 
wing, but it should be somewhere between the upper 
and the lower wing surfaces. The Ri e (rj), R te (r)), and 
Rm ( r l) are used to separate points on the upper surface 
from points on the lower surface. 

Because we know y for each grid point, we can de- 
fine a plane that passes through the grid point with a 
normal vector defined by T(rj). We must find the in- 
tersection of this plane and the curves shown in Fig. 6, 
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Fig. 6 Curves defining the backward mapping. 


T(ri) -[r - Ri e (rj)} T = 0 

(6) 

T(r))-[f-R te {r])] T = 0 

( 7 ) 

T(r,) ■[f-R m (y)] T = 0 

(8) 


Equations 6-8 must be solved for all grid points in the 
model. For a high-order NURBS curve, Eqs. 6-8 are 
nonlinear and can be solved by the Newton-Raphson 
method. The solution to Eqs. 6-8 for each ij is a set of 
three points located at the leading edge, the trailing 
edge, and the center. The %c is calculated based on 
the leading and trailing edge points. Next, we need to 
separate the grid points defining the wing model into 
upper and lower. We can connect the three points ob- 
tained from Eqs. 6-8 to form a curve that separates 
the upper surface from the lower surface. This curve 
does not have to represent the camber line accurately, 
and a wing with drooping leading edge or with highly 
cambered airfoil sections may require more than one 
R m {v) to define the curve. With this approach, it 
is possible to localize the deformation to a specific de- 
sign area by setting allowable %c m j n , %c max , r/ m j n , and 

??max • 

As the design variables (control points Pij) change, 
we can calculate the contribution from the thickness 
and camber by Eqs. 3-4. The advantage of this pro- 
cess is that the sensitivity of grid point location with 
respect to design variables is only a function of the 
B-spline basis functions, 


airfoil section incident angle. Similarly, the shear (di- 
hedral) is defined as the difference between the airfoil 
leading edge z coordinate for the root and the z co- 
ordinate at each airfoil section. If the twist angle at 
the tip is less than the twist at the root, the wing is 
said to have a washout, which could delay the stall at 
the wing tip. Also, as the wing washout increases, the 
wing load shifts from outboard to inboard. As a re- 
sult, the spa.nwise distribution of the twist, angle plays 
an important role in the wing performance. 

The SOA are used to modify the wing twist and 
shear distribution. Alan Barr presented a series of 
SOA algorithms for twisting, bending, and tapering 
an object. 13 Watt and Watt referred to these algo- 
rithms as nonlinear global deformation. 5 Sederberg 
and Greenwood extended Barr’s ideas to handle com- 
plex shapes. 17 Modified versions of these algorithms 
are presented in this paper. 

To modify the twist and shear distributions, the 
wing is embedded in a nonlinear deformation object 
referred to as a twist cylinder, that, is shown in Fig. 7. 
The twist cylinder is also used for modifying shear dis- 
tribution. The center of the cylinder is defined by a 
NURBS curve, R m (r]). The effect of deformation can 
be confined to a section of a wing by limiting the pa- 
rameter T] to vary between ? 7 m j n and r/ max . The ? 7 m j n 
can extend to the wing root, and the r/ max can go be- 
yond the wing tip. The cylinder can be twisted and 
sheared only in a plane (twist, plane) defined by a point, 
along R m (r ) ) with a normal vector of T(rj). The p\[r,) 
and p 0 (rj) are the radii of inner and outer cylinders, re- 
spectively (see Fig. 7). The deformation has no effect 
for grid points located outside of the outer cylinder, 
and the effect, of deformation is scaled linearly from 
the outer cylinder to the inner cylinder. This allows 
us to blend the deformed region with the undeformed 
region in a continuous manner. 

The angle 9{t)) defines the desired twist angle dis- 
tribution, and S(rj) defines the shearing vector. The 
9(r]) and S(rj) variables are defined by NURBS repre- 
sentation: 


OR _ OR 
idJd dPca idJd 


Njd, p (£,)Njd, q {fl) Wi d,j d 

52 Ni,p{£) 52 Nj t q(Tj)Wi t j 
i = 0 j = 0 

( 9 ) 


where id and jd are t.he indices of design variables, 
Pid,jd- Consequently the sensitivity, as seen in Eq. 9, 
is independent, of the design variables ( Pid,jd ) and the 
coordinates ( x,y,z ). Thus, we need t.o calculate the 
sensitivity with respect, t.o thickness and camber only 
at the beginning of the optimization. 


Twist and Shear 

The twist angle is defined as the difference between 
the airfoil section incident, angle at the root, and each 


%) = ~ (10) 

E N itP (n)Wi 

i = 0 

E N itP {ri)WiSi 

Sid) = ^ ( 11 ) 

E Ni,p(v)Wi 

i=0 

where 9i and Si are the twist and shear design vari- 
ables, respectively. Similar to thickness and camber 
algorithms, we use 

n = y, f(r 1 ) = (0,y,0) T (12) 
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Fig. 7 Twist definition. 


The second step for twist and shear deformation is 
to establish the forward mapping from the deforma- 
tion object (twist cylinder) coordinate system ( 77 ) to 
the model coordinate system ( x,y,z ). We use Eq. 8 
to determine 77 . Once 77 is determined, we can calcu- 
late the local 77 ( 77 ), Pi{rj), p 0 (v), T(r]), 6(rj), and S(t]). 
The point f is rotated 9{rj) degrees about R m (v) and 
sheared S. 



Fig. 8 Twist definition for a transport. 



SRtw(r]) = e(rj)p(r))[sm6(r)),0,cos6(ii)} T (13) 

SRsh(v) = e(v)S(v) ( 14 ) 

where e(rj) is a scale factor which diminishes the effect 
of deformation as we approach the outer cylinder. 


e(ri) = { 


0 

p(v)-Po 

p>(n-po 


u 


if Pin) > Poiv) 
if pi < p(rj) < p 0 (v) 
if p(rj) < p\ 


(15) 


The sensitivity of a grid point with respect to the 
twist and shear design variables is 


dR 

dOi 

dR 

dSi 


= e (v)p(v)^^-[cos9(rj),0, 


d6i 


= e(tj) 


9S(V) 

OR 


■sin 0(rf)Y( 16) 
(17) 


The term is independent of the twist design vari- 
ables Qi (see Eq. 10). However, sin $( 77 ) and cos $( 77 ) 
depend on the twist design variables and must be up- 
dated every cycle of the optimization. In contrast, the 
term d g^ is independent of shear design variables Si 
(see Eq. 11). 

Figure 8 shows the inner twist cylinder for a com- 
mercial transport. Figure 9 shows the result of twist- 
ing the wing 45° at the tip. This is a large and 
unrealistic amount of twist, but it shows the effective- 
ness of the SO A. 


Planform Parameterization 

The wing planform is typically modeled with a set of 
two-dimensional trapezoids in the x-y plane. Figure 10 


Fig. 9 Result of 45° twist on a transport. 



Fig. 10 Planform of a generic high-speed civil 
transport. 



Fig. 11 Planform definition. 

shows the planform of a generic high-speed civil trans- 
port that uses two trapezoids. As shown in Fig. 11, 
each trapezoid is defined by the root chord ( C r ), tip 
chord (C't), span ( 6 ), and sweep angle (A). From these, 
other planform parameters, such as area (A), aspect 
ratio (dR), and taper ratio (A), are defined: 

A= b -(C r + C t ), A = g- (18) 

The FFD algorithm described by Sederberg and 
Parry 14 is ideal for deforming the polygonal models. 
Like other SOA algorithms, this algorithm maintains 
the polygon connectivity, and the deformation is ap- 
plied only to the vertices of the model. The FFD 
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Fig. 12 Parallelepiped volume for free-form defor- 
mation. 



Fig. 13 NURBS volume for free-form deformation. 

process is similar to embedding the grid inside a block 
of clear, flexible plastic (deformation object) so that, 
as the plastic is deformed, the grid is deformed as well. 
Deformation of complex shapes may require several 
deformation objects. The shape of these deformation 
objects is not arbitrary. In fact, they must be three- 
dimensional parametric volumes, which could range 
from a parallelepiped as shown in Fig. 12 to a gen- 
eral NURBS volume as shown in Fig. 13. The block is 
deformed by perturbing the vertices that control the 
shape of the deformation block (e.g., corners of the 
parallelepiped). For parametric volume blocks, param- 
eters controlling the deformation are related through 
the mapping coordinates (£,»?, C)- These coordinates 
are used in both forward and backward mapping. 

Figure 12 shows a general parallelepiped defined by 
a set of control points forming three primary edges 
or directions along ?, i], and (j. The relation for a 
parallelepiped is defined as 

V, 0 = Po + hf? + n v V + »cC (19) 


grid points, are mapped to the coordinates of the 
parallelepiped, ?, ?y, and ?, as 

n n x • (r - P 0 ) 
h n x n f • (h^) 

8 < * P ° ] (20) 

x ?!<_- • (n v ) 

x h v ■ (r - Pp) 

n.£ x ii v ■ (hf) 

A grid point is inside the parallelepiped if 0 < ?, rj, ? < 
1 . 

The FFD technique based on the parallelepiped is 
very efficient and easy to implement. It is suitable for 
local and global deformation. The only drawback is 
that the use of the parallelepiped limits the topology 
of deformation. To alleviate this drawback, Seder- 
berg and Parry proposed to use nonparallelepiped ob- 
jects. 14 They also noted that the inverse mapping 
would be nonlinear and require significant computa- 
tions. 

Another popular method to define FFD is to use 
trivariate parametric volumes. Sederberg and Parry 
used a Bezier volume. 14 Coquillart at INRIA extended 
Bezier parallelepiped to nonparallelepiped cubic Bezier 
volume. 18 This idea has been further generalized 
to NURBS volume by Lamousin and Waggenspack. 19 
The NURBS blocks are defined as 


? = 
n = 
C = 


ENi, pl (?) E n J iP2 ( V ) E N kiP3 (()w itj , k P itj>k 

i = 0 j =0 k= 0 

E^*,pi(£) S Nj, p2 (v) E Nk, P 3(C) Wij'k 

i = 0 j = 0 k= 0 

( 21 ) 


where N is the B-spline basis function, and the pl- 
p3 are the degrees of N. The Pij t k are the NURBS 
control points that are related to the design variables. 
Lamousin and Waggenspack 19 used multiple blocks to 
model complex shapes. This technique has been used 
for design and optimization by Yeh and Vance 20 and 
also by Perry and Balling. 21 

The common solid elements used in FE analysis 
(Fig. 14) can be used as deformation objects. The 
mapping from the solid element coordinates is de- 
fined 22 by 


r^,vX) = Y^PiNi(i,v,0 ( 22 ) 

i 


where Po is the origin of the parallelepiped, and , n v , 
and Tiq are the unit vectors along the parallelepiped 
primary edges in ?, rj, and ? directions, respectively. 
Equation 19 defines a mapping between the deforma- 
tion object (parallelepiped) and the grid point. The 


where Ni are the FE basis functions, and Pi are the 
nodal coordinates of deformation objects, which are 
related to the design variables. The equations for the 
inverse mapping is nonlinear for all solid elements with 
the exception of tetrahedron solid elements. The solid 
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Fig. 14 FE solid elements. 



Fig. 15 Planform deformation of a transport. 


elements provide a flexible environment to deform any 
shape. Complex shapes may require the use of several 
solid elements to cover the entire domain. 

To model the planform shape, we have used hexa- 
hedron solid elements with four opposing edges par- 
allel to the 0 -coordinate. Then, the planform design 
variables are linked to the corners of the hexahedral 
elements. Figure 15 shows the initial and deformed 
model for a transport configuration. The solid lines 
represent the controlling hexahedron solid elements. 
The baseline model is on the left-hand side; and the 
deformed shape is on the right-hand side. 

As with the camber and thickness algorithms, the 
sensitivity of grid point coordinates is independent of 
the design variables ( Pidjd ) and coordinates ( x,y,z ). 
Thus, we need to calculate it only once, at the begin- 
ning of the optimization. 



Fig. 16 Implementation plan. 

Implementation 

Figure 16 shows the implementation diagram for the 
combined algorithm. The implementation starts with 
a CAD model that defines the geometry. The first two 
steps can be implemented in parallel. The first step 
is to determine the number and the locations of the 
design variables with the aid of the CAD model. In the 
second step, the grids are manually generated for all 
involved disciplines. In the third step, the mappings 
described in the previous sections are calculated for 
each grid point. In the fourth step, the new grid is 
deformed in response to the new design variables, and 
the sensitivity derivatives are computed as well. The 
third and fourth steps are completely automated. The 
first three steps are considered preprocessing steps and 
need to be done only once. 

Parameterizing Computational 
Structural Mechanics Models 

Parameterizing CFD and CSM models appear to 
be similar in nature, but the CSM model parameter- 
ization has two additional requirements. First, the 
CSM model parameterization must include not only 
the OML but also the internal structural elements 
such as spars and ribs. Second, the deformed CSM 
model must be a valid design. For example, the spars 
must stay straight during the optimization. The algo- 
rithms presented in this paper can easily handle the 
first requirement. However, if the planform design 
variables are not selected with care, the second re- 
quirement could be easily violated. To avoid creating 
invalid CSM models, planform must be parameterized 
with few hexahedron solid elements, and they must be 
aligned with major structural components such as spar 
and ribs. 
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Results and Conclusions 

The algorithms presented in this paper have been 
applied for parameterizing a simple wing, a blended 
wing body, and several high-speed civil transport con- 
figurations. Figure 17 shows the baseline and de- 
formed grids for a high-speed civil transport. The solid 
lines represent the hexahedron solid elements control- 
ling the planform variation. The parameterization 
results from this research have been successfully im- 
plemented for aerodynamic shape optimization with 
analytical sensitivity with structured 7 and unstruc- 
tured CFD grids. 9 The parameterization algorithm 
presented in this paper is easy to implement for an 
MDO application with complex configuration. The 
resulting parameterization is consistent across all dis- 
ciplines. Because the formulation is based on the SO A 
algorithms, the analytical sensitivity is also available. 
The algorithms are based on parameterizing the shape 
perturbations, thus enabling the parameterization of 
complex existing analysis models (grids). Another 
benefit of parameterizing the shape perturbation is 
that the process requires few design variables. Use of 
NURBS representation provides strong local control, 
and the smoothness can easily be controlled. 
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